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ABSTRACT 

We present the time evolution of the steUar metaUicity for SDSS galaxies, a sample 
that spans five orders of magnitude in stellar mass (10^ — 10^^ M©). Assuming the 
iBruzual fc Chariot! (|2003[ ) stellar population models, we find that more massive galax- 
ies are more metal-rich than less massive ones at all redshifts; the mass-metallicity re- 
lation is imprinted in galaxies from the epoch of formation. For galaxies with present 
stellar masses > 10^" M©, the time evolution of stellar metallicity is very weak, with 
at most 0.2 — 0.3 dex over a Hubble time- for this reason the mass-metallicity relation 
evolves little with redshift. However, for galaxies with present stellar masses < 10^" 
M0, the evolution is significant, with metallicity increasing by more than a decade 
from redshift 3 to the present. By being able to recover the metallicity history, we 
have managed to identify the origin of a recent discrepancy between the metallicity 
recovered from nebular lines and absorption lines. As expected, we show that the 
young population dominates the former while the old population the latter. We have 
investigated the dependence on the stellar models used and find that older stellar pop- 
ulation synthesis codes do not produce a clear result. Finally, we have explored the 
relationship between cluster environment and metallicity, and find a strong correlation 
in the sense that galaxies in high density regions have high metallicity. 

Key words: galaxies: evolution, galaxies: statistics, galaxies: abundances 



1 INTRODUCTION 

The quality of spectra of the observed light of unresolved 
stellar populations has reached a sufficient accuracy that 
it is now possible to make detailed studies of the physi- 
cal properties of the stellar populations in these galaxies. 
While much attention has focussed on the time evolution 
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lated quantity, the evolution of the metallicity of the star 
form ing gas in ga laxies, has not receive d as much atten- 
tion (iPanter et al.l 2003: Erb ct al. 2006; B ouche et a"l]|2007l : 
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One reason is that while it is relatively straightforward to 
measure the star formation rate from massive stars at differ- 
ent redshifts, it is not easy to determine its metallicity, which 
requires high-quality spectra. Further, in order to cover a 
large range in mass and lookback time one would need to 
observe a very large area of the sky, of the order of thousands 
of square degrees. 

One alternative approach is to use the fossil record in 



the local universe and reconstruct the star and metallicity 
histor y from this record. This however has i t s own limita- 
tions (jOcvirk et all I2OO6I : IPanter et all l2007l : iToieiro et aD 
l2007f) . Using MOPED, a tool that allows rapid extraction of 
the physical parameters of the stellar populations of galaxies 
from their spectra, our group has shown how it is possible to 
use the fossil record to obtain information abo ut the stellar 
popu lations of galaxies at large lookback times ()Panter et al.l 
l2007t) . In this work we concentrate on what can be learned 
about the metallicity of those stellar populations. 

The evolution of metallicity of the stellar populations 
of galaxies as a function of redshift is highly relevant as 
it can tell us about how the interstellar medium is being 
enriched, potentially what the initial mass function is (pro- 
vided the theoretical yields of stars are well known) and 
how this enrichment depends on environment l|Sheth et al.l 
I2OO6I) or mass, star formation history and age of the galaxy 
(' Panter et"al]|2007l ). 

Determining the metallicity history of galaxies has in- 
deed a long history. The first attempt is to use our own 
galaxy as representative and measure metallicities in its in- 
dividual stars or stellar clusters (see, for example, the review 
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bv lFreeman fc Bland-HawthornI (|2002l ')'). There are two ma- 
jor advantages: one can obtain high S/N detailed spectra, 
and the stellar populations are individually resolved. The 
main disadvantage is that our galaxy is not typical of others 
in the Universe. 

The second route is to determine the metallicity of 
galaxies from the spectrum of their integrat ed light. The 
pioneering grid of models by IWorthevI (| 19941 ) allowed this 
for simple stellar populations, i.e. elliptical galaxies. His 
models provided Lick equivalent widths in the optical re- 
gion that were sensitive to the overall metallicity of the 
stellar population. Indeed the application of this technique 
has allowed many researchers t o determine th e meta llicity 
of galaxies i n the nearby (e.g. iGallazzi et al.l (|2006l )) and 
distant (e.g. iPettinil (|200ll )') universe. The ratios of vari- 
ous emission lines can also be used to determine the metal- 
licity of the gas in a galaxy, although careful calibration 
is required llKewlev fc Dopital |2002| : iTremonti et~aLl |2004 
lErb et al.ll2006|y. 



However, Heavens et al.l l|2004l ) have shown that for an 
ensemble of galaxies it is possible to go beyond simply recov- 
ering the mass-weighted metallicity of the stars and recover 
the metallicity history in many time bi ns. This is no t neces - 
sarily true for indiv idual galaxies fe.g. lOcvirk et aP (|2006h : 
iToieiro et all l|2007l )'). which may have noisy recoveries, but 
has been shown to be unbiased when averaged over a suf- 
ficiently large sample. The recovery of the star formation 
history and metallicity for galaxies is similar to the "popu- 
lation boxes" for ne arby resolved stellar populations in local 
group galaxies (e.g. iGallart et all fioO^). 

In this paper we recover the metallicity history of SDSS 
galaxies from their spectra. We determine its evolution as a 
function of mass, how the mass-metallicity relation depends 
on the age of the stellar population and how sensitive our 
findings are to the stellar population models used. The paper 
is organized as follows: in §2 we describe the SDSS sample 
used and in §3 we describe briefly the method to determine 
metallicities. In §4 we show the derived cosmic metallicity 
history and in §5 we present a map of the stellar metallicity 
history. §6 and 7 concern the mass-metallicity relation and 
how it evolves with redshift and mass. §8 shows the impact 
of the choice of stellar population model on our findings. We 
conclude with a general discussion of our findings in §9. 



2 THE SAMPLE 



As described in lPanter et all l|2007h . the MOPED algorithm 
has been used to extract star formation and metallicity his- 
tories of a magnitude limited sample of about 300, 000 galax- 
ies drawn from the Thir d Data Release of the S loan Digital 
Sky Survey (SDSS DR3: lAbazaiian et al.l (|2005l )). 

Our main galaxy sample is determined by red apparent 
magnitude limits of 15.0 < rUr < 17.77, and we also place a 
cut on surface brightness of fj,r < 23.0. The magnitude limits 
ar e set by the SDSS targ et selection criteria, as discussed 
in lAbazaiian et al.l (|2005| ). The target criteria for surface 
brightness was Hr < 24.5, although for fir > 23.0 galaxies 
are included only in certain atmospheric conditions. In order 
to remove any bias we have therefore cut our sample at 
/ir < 23.0. At low redshifts the Sloan galaxies are subject to 
shredding - where a nearby large galaxy is split by the target 



selection algorithm into several smaller sources. To reduce 
this effect, for our star formation analysis we use a range 
of 0.005 < z < 0.34. We want to derive properties that 
avoid dependence on the SFH of the galaxies, and hence 
wish to avoid inverse Vmax weighting where possible. We 
therefore choose mass cuts around a narrow redshift range 
of 0.1 for the deriving the cosmic metallicity evolution, and 
show colors normalised to the total number of galaxies in 
that mass range in the mass-metallicity diagrams. 



3 DETERMINING THE METALLICITIES OF 
GALAXIES 

In order to determine the optimal multi-population fit for 
each spectrum it is necessary to allow as free a star for- 
mation history as possible. A parametrization of 11 SSPs 
spaced logarithmically in lookback time, each with inde- 
pendent and variable metallicity, was allowe d along with a 
one p arameter slab extinction based on the ICordon et al.l 
l|2003t ) LMC curve. In order to fully assess the resulting 23 
dimensional likelihood surface it was necessary to u se the 
MOPErFl ltHeavens et al, ,200(]|: iReichardt et allbOOll ) algo- 
rithm co upled with a Mar kov Chain Monte Carlo search 
method jPanter et al.l 2003h. develope d to use high spectral 
resolution models ( Panter et al.|[2007l ). For each galaxy the 
best solutions from a series of 100 randomly-seeded conju- 
gate gradient searches were used to seed a Markov chain of 
10^ steps. The individual search traces were examined to 
confirm adequate exploration of the parameter space and 
estimate solution convergence. 

In our approach, for each galaxy we describe its star 
formation by 11 lookback time bins (spaced logarithmically) 
and for each of these we determine the mass fraction trans- 
formed into stars and its metallicity. The first interesting re- 
sult, before looking at the metallicity history of galaxies, is 
to look at the overall metallicity of the stars in SDSS galax- 
ies today. It is important to realize that there must be some 
weighting of the metallicity values recovered - the metallic- 
ity of a population that contributes very little light to the 
final spectrum is not constrained, so should be discounted. 
Since the MOPED method recovers the complete star for- 
mation history of a galaxy it is possible to calculate the 
average metallicity with different weights. Traditional lick 
index analysis considers the light-weighted metallicity, but 
with knowledge of the SFH and IMF it is possible to deter- 
mine present mass fraction (PMF) or original star formation 
fraction (SFPfl average metallicities weighted by either. 

In Fig. [T] we illustrate how a synthetic galaxy spectrum 
responds to changes in metallicity history, by comparing a 
synthetic spectrum created using the mean star and metal- 
licity history of SDSS galaxies recovered by MOPED with 
models that have the same star formation history but differ- 
ent metallicity histories. It is clear that changing the metal- 
licity has a significant effect on both the continuum level 
and depth of absorption features. 



^ Massively Optimized Parameter Estimation and Data Com- 
pression 

^ i.e., the fraction of total star formation over the galaxy's history 
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Figure 1. A synthetic spectrum corresponding to the mean star 
formation history of SDSS galaxies, and the effect of changing 
the metalhcity history while maintaining the same star forma- 
tion history. The lower panels show the difference between the 
continuum-subtracted mean spectrum and a model produced us- 
ing the different continuum-subtracted alternate metallicity mod- 
els, illustrating that changes in metallicity history can be seen in 
both the low and high frequency spectral signal, fn this figure, 
increasing from solar refers to an initial metallicity of Zq in the 
oldest bin rising to 2.5Zq in the youngest and decreasing from 
solar an initial metallicity of Zq falling to O.O2Z0. 



We can distinguish regions in wavelength that are driv- 
ing the differences coinciding with H, Ca, Fe and Mg absorp- 
tion lines. The continuum shape blueward of 4000A also has 
an impact, but in this plot we are showing only the high 
frequency signal of the residuals normalized to the contin- 
uum. While traditional indices based approaches will con- 
sider some subset of this data, by analyzing the full spec- 
trum we are able to determine the metallicity of individual 
galaxies more accurately. 



4 THE COSMIC METALLICITY HISTORY 

In order to combine the metallicity histories of many galaxies 
at different redshifts it is necessary to estimate metallicity 
histories in a common time frame. In previous work we have 
formed complex transform functions to redistribute a galax- 
ies SF in a new set of bins, but in this case we instead chose 
a narrow redshift range of galaxies — 0.1| < 0.01). For our 
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Figure 2. Bulk metallicity of stars in SDSS galaxies for different 
stellar masses. This is the average metallicity of all stars present 
at a particular redshift, regardless of age. To avoid possible errors 
from shifting bins we have chosen a narrow redshift range centred 
at z = 0.1. The number of galaxies in each mass range is given 
in the inset. Note that more massive galaxies are more metal 
rich than less massive galaxies. The enrichment is very fiat as a 
function of redshift for the more massive galaxies, as the bulk of 
their stars formed early and the tiny amount of recent SF has not 
been sufficient to change the bulk metallicity. The dotted lines 
correspond to a Icr bootstrap error. 



common bins we use the original time bins (see lPanter et al] 
liooj)) with i = at 2 = 0.1. This gives a sample of 42160 
galaxies. 

The simplest interpretation of the metallicity history is 
that of the bulk of stars in a galaxy at a given time. To 
calculate this for each bin we weight the metallicity in the 
bins up to that point by the mass of stars remaining from 
that bin (ie, we include recycling). This traces the average 
metallicity of stars in a galaxy over its evolution, and is 
shown in Fig. (2] 

The individual metallicity histories of each galaxy can 
also be combined to give the mass-weighted cosmic metal- 
licity history of the star forming gas. To calculate this we 
average the metallicities for a given bin, weighted by the 
mass of that population. We must also impose a cut to en- 
sure that the metallicity of populations which have a very 
small contribution to the spectrum are not considered, as 
these have poorly determined metallicity and give results 
dependent on the metallicity prior. To avoid this contam- 
ination we add a further cut: to contribute to the average 
of a particular time bin a galaxy must have > 25% of its 
spectral flux in that bin (see Fig. 

Weighting by bin mass obviously biases the results to- 
wards high-mass galaxies. We also present results using a 
weighting scheme that uses the fraction of the total spectral 
flux. Weighting by mass means that there is a slight bias 
towards older populations, which may not be as well deter- 
mined as for a given mass as the younger populations give 
much more light. To calculate the cosmic metallicity history 
the metallicities in the bins are then averaged, weighted by 
the fractional contribution of the bin to the galaxy luminos- 
ity. In this weighting scheme a galaxy's mass is irrelevant. 
The corresponding cosmic metallicity history is shown in 
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Figure 3. Average metallicity for different stellar masses, as in 
Fig. 1 except for star forming gas rather than stars. 
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Figure 4. Same as Fig. |3] but for a total fractional spectral flux 
weighted average. 



Fig. |4] and is very similar to the flux-weighted history (Fig. 

[31). 

The flrst thing we learn is that most metals are locked 
in the most massive galaxies and that the metallicity evo- 
lution for each mass range is relatively flat as a function of 
redshift, with variations of 0.1-0.2 dex for each mass range. 
Note that a similar trend is found for a light-weighted plot. 
Our results compare favourably with the trends observed at 
high redshift by Maiolino et al. (2007). If we look at their 
Fig. 3 and compare their abundance as a function of red- 
shift for a stellar mass of 1 x 10^" Mq with our green line 
in Fig. 1, 2 and 3 we find that their abundance values at 
z = 3,0.07 and are log(Z/Zo) = -0.54,-0.1,0.1 respec- 
tively, assuming that 12-|-log(O/i?)0 = 8.66. This compares 
with our derived values, for the same redshift intervals, of 
—0.5, —0.2, 0.0, which is within 0.1 dex of their values, which 
is remarkable. Any conclusions are subject to model choice, 
as we explore in §8. 



Table 1 . Pearson correlations of metallicity and number of cluster 
galaxies in a cell 



z 


Smoothing 


Radius(°) 







2 


3 


4 


0.187 


0.091 


0.118 


0.104 


0.079 


0.456 


0.101 


0.167 


0.171 


0.171 


1.21 


0.131 


0.213 


0.212 


0.189 


3.64 


0.127 


0.144 


0.136 


0.112 



5 MAPPING THE METALLICITY 
EVOLUTION OF THE UNIVERSE 

The SDSS-DR3 spectroscopic footprint covers 3732 sq. deg., 
roughly 10% of the sky. We use the metallicity history of the 
galaxies to create maps over this area of the enrichment his- 
tory at different epochs. We use the HEALPi^lfl algorithm to 
determine equal area patches on the sky, and calculate the 
mass weighted average metallicities for each patch and time 
bin as before. Figure [S] shows the mass-weighted metallic- 
ity maps for our four highest redshift bins, smoothed with 
a boxcar filter of radius 2°. Over-plotted are the locations 
of t he Brightest Clust er Galaxies from the SDSS C4 cata- 
log l|Miller et al.ll2005l ). used to represent the distribution of 
cluster galaxies on the sky. It is clear by eye that in many 
regions the crosses follow the regions of higher metallicity 

For areas of the footprint where cluster galaxies ex- 
ist, a cross-correlation analysis between mass weighted av- 
erage gas metallicity and number of cluster members in 
cells reveals strong correlation between the three oldest bins 
z = 0.456, 1.21, 3.64 and the number of cluster galaxiefl 
(see Tables [ 1] and [H). It would appear that, as stated in 
ISheth et al.l ( 20061 ) . metallicity is strongly correlated with 
environment - this can be interpreted as the seeds of clusters 
being the seeds of metal enrichment in the universe. Note 
that the overall level of enrichment in the map ai z = 3.64 
is very homogeneous at around the solar metallicity value, 
while ai z — 0.187 there is much more variation. If we as- 
sume that metallicity is an indicator of environment, this 
offers a tantalisingly glimpse of the growth in the influence 
of dark matter structure, only visible by examining the huge 
volume at high (temporal rather than spatial) redshift of- 
fered by the fossil record. 

It will be interesting to cross-correlate these maps with 
the upcoming SZ experiments which will detect clusters of 
galaxies at higher redshift. The larger variations in the low- 
est redshift map probably reflect the change in sampling to 
less massive galaxies, as these are those in the sample which 
are likely to have a high fraction of younger star formation. 



package are available from 



3 Details of the HEALPix 
http : //healpix . j pi . nasa . gov 

* NB: We correlate with cluster members, not BCGs 
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Figure 5. HEALPix projection of the mass weighted gas metaUicity for several look back bins averaged over galaxies in the redshift 
range 0.0 < z < 0.1. Galaxies only contribute to a bin if > 25% of their star formation occurs in that bin. The overlaid crosses correspond 
to the bright cluster members in the C4 catalog. The metallicities are boxcar smoothed over a 2° radius. The size of the smoothing patch 
is shown in the bottom right of each panel 
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Table 2. Spearman Correlations of metallicity and number of 
cluster galaxies in a cell 
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Smoothing 


Radius(°) 







2 


3 


4 


0.187 


0.077 


0.097 


0.078 


0.061 


0.456 


0.084 


0.179 


0.182 


0.167 


1.21 


0.127 


0.274 


0.273 


0.246 


3.64 


0.169 


0.201 


0.189 


0.154 



6 THE AVERAGE MASS-METALLICITY 
RELATION 

We now turn our attention to the local mass-metallicity rela- 
tion for the SDSS galaxy population to understand its origin 
and its time evolution. 

Fig. |6] shows the mass-metallicity relation for 312, 815 
galaxies in the DR3 Main Galaxy Sample. In this case we 
calculate the mass fraction weighted metallicity, ie for an 
individual galaxy the metallicity is calculated by weight- 
ing the bin metallicities by the observed fraction of mass in 
that bin. The three solid lines represent the 16th, 50th and 
84th percentiles of the distribution. The first thing to note is 
that there is a clear mass-metallicity relation: more massive 
galaxies harbour stars with higher metallicity. The average 
metallicity of an L, galaxy is solar. For lower masses, the 
metallicity decreases approximately by 0.5dex for every dex 
in mass. 

For galaxies with stellar masses of about 10^ Mq, the 
average metallicity of the stars is 0.1 the solar value. Note 
that for masses larger than 10^^ Mq there is a flattening of 
the mass-metallicity relation. The maximum value we obtain 
for the mass-metallicity relation is I.IZq. The spread in the 
relation is also smaller at higher masses (0.15 dex) and grows 
at smaller masses (0.5 dex). 

There is a break around M, — 10^° M©, below which the 
dispersion around the median value increases, and a much 
wider range of metallicities is recovered for a given mass. 
This can be interpreted as consistent with the turnover p oint 
noted in iKauffmann et al. 1 l|2003l ): Ijimenez et all l|2005h . 

We can obtain a reasonable fit to the median line using 
a tanh function over the mass range where cells contain more 
than 2000 galaxies (8.8 < logM. < 11.8) of the form 



log — — = A + B tanh 

Zq 



log M, - log Mc 



(1) 



where logAf^ = 9.66, A = 1.04, A = -0.452, B = 0.572 
and masses are in solar masses. This formula is given purely 
for convenience, and should not be taken to reflect any un- 
derlying physical motivation or be extrapolated outside of 
the fitted mass range. It is of interest though that the break 
at log Mc ~ 9.66 neatly corresponds with the increase in 
dispersion noted above. 

We can now compa r e this relation with previous work. 
Recently iGallazzi et al.l (j2005l ) , using absorption indices in 
the spectra of SDSS galaxies, have determined the local 
mass-metallicity relation. Their approach is different from 
ours since they concentrate on specific absorption features 
of the spectrum, and their results are weighted towards the 
most luminous populations. For most old galaxies this will 
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Figure 6. Mean metallicity mass relation using the present mass 
fractions as weights. The solid green lines represent the 16th, 50th 
and 84th percentiles of the distribution and the light blue a tanh 
fit to the results . The dashed and dot -dash lines a re the same per- 
centiles for the IGallazzi et al.l ((2005) (red) and iTremonti et al.l 
|2004|) (orange) results. The distributions are normalized over 
each column of mass, the underlying galaxy mass distribution 
is shown in the lower histogram 



be very similar to our own result, but because we know the 
contribution to the spectra of each component we can recal- 
culate the average metallicity of galaxies with similar weight- 
ing (Fig. [7]). The agreement is very good over the match- 
ing mass range, although seems to be breaking down at the 
low-mass end. This is not surprising because due to down- 
sizing of the galaxy population, it is the low-mass galaxies 
that dominate star formation today. This of course is due 
to the fact that for younger galaxies with ongoing star for- 
mation it is the younger stars that are more luminous, even 
though by mass the galaxies are almost certainly dominated 
by older stars (Panter et al. 2007). The otherwise excellent 
agreement reinforces the fact that we are not introducing 
degeneracies in the overall quantities by recovering the star 
and metallicity history in 11 time bins. In Fig. [7] we show 
the same plot as in Fig. [6] but this time luminosity- weighted 
(strictly, by spectral flux). Note that in this case there is 
again good agreement, slightly better than before for lower 
masses. 

We have al so overplotted t h e mas s-metallicity relation 
determined by iTremonti et al.l (|2004[ ) from emission line 
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Figure 7. Mean metallicity mass relation using the luminous 
fractions as weights. The solid line represents the 16th, 50th and 
84th percentiles of the distribution. The dashed and do t-dash 
lines are the same percentiles for the Galla zzi et al 

] \200± (red) 

and lTremonti et al.l | |2004| ) (orange) results. The distributions are 
normalized over each column of mass, the underlying galaxy mass 
distribution is shown in the lower histogram 

measurements, and a clear offset can seen. In the next sec- 
tion we take advantage of MOPED ability to recover the 
time evolution of the metallicity and star formation history 
to explain the difference. 



7 RECOVERING THE MASS-METALLICITY 
FOR YOUNGER POPULATIONS 

Since the MOPED algorithm recovers both the star forma- 
tion and metallicity history of a galaxy it can be used to in- 
vestigate the evolution of the metallicity for a given galaxy 
- a task impossible with simple indices. In particular we can 
look at only galaxies with recent star formation and exam- 
ine how the metallicity of the star forming gas has evolved 
from that which produced the majority of the stars at earlier 
ages. 

By using only the younger populations (< IGyr) of 
galaxies for the analysis we can reconstruct the informa- 
ti on that would be prob ed by the emission line analysis (as 
in iTremonti et al.l (|2004 l). Note that our analysis explicitly 
removes strong emission lines - all our information comes 
from continuum shape and absorption features. Excluding 




Present Stellar Mass / [M^] 

Figure 8. Intermediate-age mass-metallicity relation. To com- 
pose this plot we have considered only those galaxies with more 
than half their spectral flux in populations < IGyr. The metal- 
licity is the average of the younger bins weighted by their frac- 
tional spectral flux. The solid line represents the 16th, 50th and 
84th percentiles of the distribution. The dashed and do t-dash 
line s are the same p ercentiles for the ICallazzi et al.l (|2005| ) (red) 
and lTremonti et al.l 1 12004 ) (orange) results. The distributions are 
normalized over each column of mass, the underlying galaxy mass 
distribution is shown in the lower histogram 

AGN activity, the emission lines reflect the metallicity of 
nebular gas around star-forming regions that has been ex- 
cited by the UV emission of young stars. Fig. [8] shows the 
mass-metallicity relation of the star-forming gas as derived 
indirectly from our analysis. It is important to use only 
metallicities which are well determined - for this analysis we 
use the 27915 galaxies for which over half the spectral hght 
comes from these younger populations, and plot the metal- 
licity weighted by their fra ctional spectral f lux. N ote that 
now the agreement with the lXremonti et aP (|2004h result is 
much better and our analysis shows that their metallicities 
are dominated by an intermediate ~ 1 Gyr stellar popula- 
tion in sub— L* galaxies. In Fig. |9] we extend this analysis 
by comparing directly the metallicities determined from ab- 
sorption and emission of galaxies which have more than half 
of their flux in populations < IGyr. The figure shows good 
agreement, confirming that given sufficient flux in those pop- 
ulations the MOPED algorithm using absorption features 
can determine accurate metallicities. The offset could indi- 
cate a systematic difference in metallicity between the birth 
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Emission Metollicily / 1 2 + log(0/H)-B.69 

Figure 9. A comparison of metallicities determined by absorp- 
tion and emission analysis. This plot considers only those 14038 
galaxies with more than half their spectral flux i n populations < 
IGyr that also have a metallicity reported in the lTremonti et al.l 
1I2OO4I ) sample. The absorption metallicity is the average of the 
younger bins weighted by their fra ctional spectral f lux an d the 
emission metallicity is taken from iTremonti et al.l (|20o3) ■ The 
blue diamonds show the median value of absorption metallicity for 
a given emission metallicity, the black line is a linear fit through 
all points and the red dotted line corresponds to a 1:1 relation. 
See text for a discussion of the offset. 



cloud (responsible for the emission features) and the young 
stars, or possibly a calibration offset in the models. 

We can extend this analysis to galaxies younger than 
the 1 Gyr bin, although the number of galaxies with suffi- 
cient star formation in these youngest age bins is very small 
and determination of metallicity in these younger stars with 
weaker absorption features is more difficult than in older, 
cooler stars. Fig. [10] shows the mass- metallicity relation of 
these younger populations and it is clear that any relation- 
ship has broken down. First, the stellar mass range is much 
smaller than before because of downsizing in the galaxy 
population. Only galaxies with stellar masses below 10^*^ 
Mq contribute significantly to the young population metal- 
licity and the mass range is only of 1.5 dex. As pointed 
out earlier the mass-metallicity relation is fiat (at least for 
M < lO^'^'^M©), with no dependence on mass and an aver- 
age metallicity of 0.2Z/Zq . Why galaxies below this mass all 
have the same metallicity is something we will investigate in 
future work, but one could envision environment as a major 
driver, where the gas in filaments that feed the galaxy has 
been enriched to a certain level. 



8 MODEL DEPENDENCE OF THE 
METALLICITY HISTORY 

The results shown so far have all been generated us- 
ing the current s tanda rd in stellar modelling codes 



iBruzual fc CharlotI \200^ ). It is important however to con- 
sider the effect of changing the models used; we make 
no statement as to which set is "right", and emphasize 
that the set of gal axies chosen for comparison, detailed in 
iPanter et all (|2007l '). is relatively small (808). 

We have investigated the impact of model choices 




- 300 S 



-200 



9.0 9.5 10.0 10.5 
Present Stellar Mass / [M^] 



Figure 10. Youngest age mass-metallicity relation. To compose 
these plots we have considered only those galaxies with more than 
half their spectral light in populations < IGyr. The metallicity 
is the average of the younger bins weighted by their fractional 
spectral flux. The solid line represents the 16th, 50th and 84th 
percentiles of the distribution. The dashed and dot -dash lines 
are the same p ercentiles for the ICallazzi et al.l ||2005| ) (red) and 
iTremonti et al.l l i2004| ) (orange) results. The distributions are nor- 
malized over each column of mass, the underlying galaxy mass 
distribution is shown in the lower histogram 



by recovering the metallicity history with a suite of dif- 
feren t stellar populatio n models: SPEED (Timenez et al.l 



, PEGASE (Fioc fc_Rocca-Volmerangelll99^) , BC9 3 



(Br uzual fc CharloTTiggal'). Marasto n (|Marastonl 120051 ') . 
GALAXEV iBruzual fc Charlotll2003l ). CB07 (As BC03 but 
an with improved treatment of TP-AGB stars). Our refer- 
ence model is the BC03 at a spectral resolution of SA. 

In Fig. [Tl]we show the differences in average metallicity 
of stars as a function of redshift. It is clear that our conclu- 
sions regarding higher-mass galaxies having higher metallic- 
ity is robust whichever model is used. The deviation shown 
by the SPEED mo dels refiects the fact that, as shown in 
jPanter et alll|2007l ). the star formation history recovered by 
this model is broader, with overall SFR peaking at z < 1. It 
should also be noted that the SPEED models allow a wider 
range of metallicity values than other synthesis codes. 

In Fig. 1121 there is apparently some issue at z=0.2 (f ~ 1 
Gyr), where almost all the models give very similar metal- 
licity. With a population of galaxies taken at 2 = 0.1, this 
corresponds to a population of stars of age 1 Gyr, in some 
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Figure 11. Cumulative metallicity of stars in SDSS galaxies re- 
covered using different stellar population models. The lower panel 
illustrates the difference between low and high mass galaxies. 



Figure 12. Metallicity of star forming gas in SDSS galaxies re- 
covered using different stellar population models. We caution that 
the region lower than z = 0.3 (dotted) has relatively little flux, 
and is not used for our earlier analysis. 



respects a difflcult population to identify as no part of the 
spectrum is domi nated by this population for usual star for- 
mation histories franter et "al]l2007l : [Mathis et al.ll2006t) . As 
this age reflects the point for the average galaxy where flux 
per population is at a minimum, we caution against over- 
interpretation. Although it may be a systematic failing of 
the libraries used to calibrate the various sets of models it 
seems more likely that what is being reflected is the prior 
introduced by a fiducial metallicity. Alternatively, it could 
be that a more sophisticated treatment of the post-AGB 
stages, as included in the more modern models, removes the 
problem. Whichever set of models is chosen, the redshift 
evolution of the stellar metallicity in the SDSS population 
is fairly flat, with a variation of about 0.5 dex. 

Of particular interest are the differences between the 
Maraston, SPEED and CB07 models. These three each in- 
clude a treatment of mass loss along the thermal-pulsating 
AGB phase of stellar evolution. The SPEED mo dels incor- 
porate this effect using an empirical prescriptio n l|jorgensenl 
Il99ll : IJimenez et aLlll998l : I Jimenez et aLll2004l ). If we focus 
on these models which include mass loss (Maraston, SPEED 
and CB07) then the discrepancy in recovery of metallicity 
is of about 0.1 dex beyond z=l, but diverges by 0.5 dex 
at lower redshifts. Clearly the metallicity is more sensitive 
to the mode l choice than the star formation history (see 
IPanter et all \2m% ) and the error budget in determining 



the metallicity history is dominated by the choice of model. 
Further work is needed to calibrate the models with higher 
accuracy. 

Figure [T3] shows in different panels the effect of vari- 
ous systematics. It is clear from Fig. I13f a) that changing 
the photometric reduction pipeline will have an effect on the 
metallicity values recovered (the SPEED model was used for 
this comparison, GALAXEV for the following). The changes 
from DRl to DR3 are clear, but it appears that PCA skyline 
(|Wild fc Hewettll2005l ) cleaning does nothing to change the 
recovered metallicity history. Fig. I13f bl makes it clear that, 
for a model calculated and calibrated at SA but rebinned to 
20A it is still possible to recover the same metallicity his- 
tory. It is interesting however to note that this is the only 
change which has any effect on the 1 Gyr bump evident in 
most of other stellar population models, but further analysis 
is required (with many more galaxies and rebinning resolu- 
tions) to resolve this issue. Figures llSf cl and llST d) show 
that changing IMF and dust extinction curve do not alter 
the metallicities recovered by a large amount. The appar- 
ently large variations in the right hand panels can be ex- 
plained by the changes in galaxy mass, caused by changing 
either the IMF or extinction curve, moving galaxies between 
the two samples. 
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Figure 13. Effect of several systematic effects on the determination of the metallicity from the spectra of the stellar populations of 
galajdes (left), and how they effect conclusions based on mass (right). The different panels show the effect of different systematics: Fig. 
(a) shows the effect of different pipelines in the data reduction of the SDSS spectra, the error is comparable to the random error, although 
it is clear that there is a difference between the DRl and DR3 reductions. Fig. (b) shows the effect of spectral resolution in the models; 
it would not be surprising to see a difference when the spectral resolution is increased, but this is not evident. This could be due to the 
fact that the underlying model (BC03) is calibrated at SA and rebinned to 20A. Fig. (c) shows the difference by choice of IMF - the 
overall change in recovered metallicity of the star forming gas is small, but the large changes in mass around the cut (due to the strong 
dependence of mass on the IMF) give a misleading effect on the right hand plot. Fig. (d) shows the dependence on dust model choice, 
and again the variation on the right hand figure is due to mass differences rather than a large difference in recovered metallicity. 
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9 CONCLUSIONS 

In this paper we have presented the cosmic evolution of the 
bulk metallicity of galaxies, demonstrating a clear trend of 
downsizing which is robust to the selection of stellar popu- 
lation model to determine the metallicity of the galaxy. We 
have extended this to explore the evolution of metallicity 
of star forming gas, which also follows a downsizing sce- 
nario but is more dependent on the choice of stellar model. 
The trend shows how massive galaxies have higher metallic- 
ities at earlier times while less massive galaxies have lower 
metallicities at earlier times but reach the same metallic- 
ity at about 1 Gyr from the current time. We have shown 
a clear correlation between metallicity and cluster environ- 
ment, and will perform more correlations as further datasets 
become available. 

We have reconciled the disparity between the mass- 
metallicity relations recovered from nebular emission line 
spectra and stellar absorption methods. As one might ex- 
pect, the absorption method is entirely appropriate for esti- 
mating the metallicity of older, established galaxies, while it 
may fail for galaxies which are undergoing recent star forma- 
tion where the most luminous population is almost certainly 
not the most massive in stellar mass. We have verified that 
the fossil analysis technique, which here excludes emission 
lines, can be used simultaneously to uncover the metallici- 
ties of both young and old populations by producing a mass- 
metallicity relation that accurately maps the results given 
by both emission and absorption line diagnostics. 

It is clear from the comparison of stellar models that 
the metallicity history determined for the star forming gas 
is highly model dependent. Although we use models with a 
wide range of publication dates, it is apparent that the three 
most modern show the largest deviation. We acknowledge 
that alpha enhancement, not included in any of the stellar 
population models we consider, could affect the metallicities 
we recover. 
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